Exchange interaction of magnetic impurities in graphene 
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The indirect exchange interaction between magnetic impurities localized in a graphene plane is 
considered theoretically, with the influence of intrinsic spin-orbit interaction taken into account. 
Such an interaction gives rise to an energy gap at the Fermi level, which makes the usual RKKY 
model not applicable. The results show that the effective indirect exchange interaction is described 
by a range function which decays exponentially with the distance between magnetic moments. The 
interaction is also shown to depend on whether the two localized moments belong to the same 
sublattice or are located in two different sublattices. 
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I. INTRODUCTION 

Physical properties of two-dimensional (2D) graphene 
planes are a subject of current interest due to their unique 
transport properties Some of these properties orig- 
inate from peculiarities of the corresponding electronic 
structure, or more specifically from the peculiar Fermi 
surface of a 2D graphene plane, which consists of six 
points located at the corners of the corresponding hexag- 
onal 2D reciprocal latticeA^ Only two of these points, 
however, are nonequivalent. Apart from this, the low 
energy excitations at these points can be approximated 
by the energy dispersion relation linear in the wavevec- 
tor. Such electron states can be described by the Dirac 
equation (therefore, the corresponding points of the Bril- 
louin zone are also termed as the Dirac points). All these 
features of the electronic structure make graphene very 
attractive from both theoretical and experimental points 
of view. This is because graphene is an ideal natural sys- 
tem for testing various theoretical models of 2D electron 
transport, including for instance the quantum Hall effect, 
or the effects due to the geometrical Berry phasei^ 

As concerns magnetic properties of graphene, these 
are still a subject of discussion in the relevant litera- 
ture. Some experiments indicate on the existence of a 
spontaneous magnetic moment of the graphene planes. 7 
However, physical origin of the spontaneous magnetiza- 
tion is not well understood. Several distinctly different 
mechanisms leading to ferromagnetism of graphene have 
been proposed, including the instability of the paramag- 
netic phase due to electron-electron interaction, 8 indirect 
RKKY exchange coupling between magnetic moments of 
structural defects via mobile electrons^ and others^ 

The indirect exchange interaction between local mag- 
netic moments is generally determined by electron ex- 
citations near the Fermi energy. Assuming that the 



Fermi momentum in a graphene plane is zero and the 
excitations are gapless, a non-oscillatory indirect ex- 
change interaction of the RKKY type has been found in 
Ref. |9(. Moreover, the corresponding exchange integral 
was found to be ferromagnetic, with the range function 
decaying as i? -3 with the distance R between magnetic 
moments. However, the intrinsic spin-orbit (SO) interac- 
tion in graphene has been neglected in Ref. |9j. Such an 
interaction in graphene opens an energy gap at the Fermi 
levels and makes the usual RKKY mechanism not rel- 
evant due to the absence of electrons at the Fermi level. 
Therefore, another mechanism of the indirect interaction 
between magnetic moments of structural defects should 
be developed. 

In this paper we reconsider the indirect exchange in- 
teraction in graphene, taking into account the role of in- 
trinsic SO interaction. If the SO gap is relatively small, 
the excitations across the gap contribute to the indirect 
exchange interactioniii* 1 ^ We show that the exchange 
coupling is ferromagnetic, and the corresponding range 
function decays exponentially with the distance between 
magnetic moments - provided the chemical potential lies 
in the gap. Moreover, the spin-spin coupling is shown 
to be described approximately by the isotropic three- 
dimensional (3D) Heisenberg Hamiltonian. Apart from 
this, the interaction is shown to be dependent on the 
location of magnetic impurities - the interaction for im- 
purities located in the same sublattice is different from 
that between impurities located in different sublattices. 

II. MODEL AND EXCHANGE INTERACTION 

Indirect exchange interaction between magnetic mo- 
ments of structural defects (for instance impurities) is 
determined by the electronics structure of the relevant 
system. To describe electron states in a graphene plane, 
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we assume the Hamiltonian of non-interacting electrons 
in the following iorm—: 



H = fiv [r x p z k x + Tyk y ) + Ar z p z a z 



(1) 



where three types of the Pauli matrices; r, p, and <r, 
have been introduced. These matrices operate in differ- 
ent spaces, so their products in Eq.(l) should be under- 
stood as the direct matrix product. The matrix r acts 
in the (real) space of two nonequivalent sublattices. The 
matrix p, in turn, acts in the (reciprocal) space of two 
nonequivalent Dirac points, whereas the matrix a oper- 
ates in the electron spin space. Finally, the constants v 
and A in Eq.(l) are the velocity of electrons near the 
Fermi level (in the limit of no intrinsic SO interaction), 
and the parameter of SO interaction, respectively. 

As we have already mentioned above, the Hamiltonian 
(1) operates in the 8-dimensional space, and describes 
electrons in the vicinity of two nonequivalent Dirac points 
located at the corners of the Brillouin zone, ko and — ko, 
with k = (27r/3a , 0) and a denoting the lattice con- 
stant. The wavevector k in Eq. 1 is measured from the 
location of a given Dirac point. 

The energy spectrum of the Hamiltonian (1) consists 
of four-fold degenerate dispersion curves given by 



e fe = ±(A 2 



ftW) 1/2 . 



(2) 



As follows from Eq. (2), the SO interaction opens a gap 
in the electronic states around the Fermi level, which 
extends from — A to A. 

We assume that the interaction of 2D electrons with 
magnetic impurities located at Rj (i = 1,2,...) in the 
graphene plane can be described by the following Hamil- 
tonian: 



H t = [1 + (-i)'T,] fa -a) 6 



(3) 



where Si is the spin moment of i-th impurity, go is the pa- 
rameter of exchange coupling between the localized spin 
and band electrons, and n is the areal density of the host 
atoms in the graphene plane. The parameter s — 0, 1 dis- 
tinguishes the two (A and B) sublattices - if a given im- 
purity belongs to the sublattice A (B) then s — (s = 1). 
This form of interaction accounts for the exchange cou- 
pling that is sensitive to the location of the magnetic im- 
purities in the two sublattices. However, the interaction 
(3) does not take into account the intervalley transitions. 
We consider first this simplified magnetic interaction, as- 
suming that the rate of intervalley transitions is small. 
Later on we will briefly analyze the effect of intervalley 
scattering on the indirect exchange interaction. 

The Green function corresponding to the Hamiltonian 
(1) has the form 

r / t_ \ _ e + hv (r x p z k x + Tyky) + Ar z p z a z 

The corresponding Green function in the energy- 
coordinate representation can be derived from the above 



equation by integrating G e (k) over the momentum in the 
vicinity of each valley, 



G e (±R) = |^ e ±^o+k).R G£(k); 



(5) 



where t = ±1 corresponds to the two valleys located at 
±ko, respectively. 

Taking into account Eqs. (4), (5) and following Ref. 0] 
one obtains 



G e (±R) = e 



±itkn-R 



i (e + AT z p z a z ) 
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where H„ (z) are the Hankel functions which decay ex- 
ponentially in the upper half-plane of the complex argu- 
ment zi^ The Green function defined by Eq. (5) is the 
analytical function in the complex plane of e, except for 
two cuts along the real axis from — oo to —A and from 
A to +oo. 

In a general case, the indirect exchange interaction is 
mediated by inter- and intra- valley electron excitations. 
As already mentioned above, we consider first the intra- 
valley contribution and then discuss the role of inter- 
valley transitions. In the loop approximationMii& we find 
the interaction energy of two magnetic impurities Si and 
S2 separated by a distance R in the form 



Eint(R) — Wap{R) Si a S2/3, 



(7) 



where 



1 + (-1) s t z apG s (-K). (8) 



Making use of Eqs. (5) and (7) and calculating the trace, 
we find 
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where the integral is along the real axis above the cut. 
After rotating the integration contour to the complex e- 
plane along the imaginary axis, as explained in Ref. |14| . 
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we find from Eq. (9) 

- aV / 



7r 3 n 2 ?i 4 i> 4 Jo 



d£{5 8 



I ( fl ^ +A2 ) [(? + A 2 ) 6 a p - A 2 6 az 5 l 3 z ] 

+<w (i - s ss ,) (e + a 2 ) k\ ( R y/e+^ 1 



,(10) 



where K v {z) are the MacDonald functions related to 
Hi X) {z) byi£ 

Kv (z) = — e^HP (ze™/ 2 ) . (11) 

When the spin-orbit interaction is taken into account, 
A ^ 0, the range function given by Eq. (10) scales with 
the characteristic length Rq = hv/A. The derived ex- 
change interaction is generally anisotropic. The interac- 
tion between in-plane spin components (i = x, y) is given 
as 



< (R) 

where 



■^{^ + (1-^)1^)}, 



(12) 



/ 0) iOR) = J™dx(x 2 + l) tf^*^±±j. (13) 

On the other hand, the interaction of the out-of-plane 
spin components is expressed via the formula 



with 



(14) 



(15) 



The interaction described by Eqs. (12) and (14) is for- 
mally valid under the condition R> a, with a being the 
lattice constant. 

The asymptotics of MacDo nald functions at 
large argument^ K v (z) ~ y / 7r72ze~ 2: , leads to 
the exponential decrease of interactions, Iq^R) ~ 
R- 3 / 2 exp{-2R/R ) and Z 2 (-R) ~ R~ 5/2 exp(-2i?/i? ) 
for R ^> Rq. Thus, the value of i?o determines an 
effective range of the exchange interaction. 

The magnitude of the indirect exchange interaction de- 
pends on the relative position of the spins in the graphene 
sublattices. As mentioned above, the exchange interac- 
tion is described by an anisotropic Heisenberg model. 
However, Iq(R) ~ ^(-R) , so the anisotropy is weak 
and one may approximate the exchange interaction by a 
3D-isotropic Heisenberg Hamiltonian, irrespective which 
sublattice the magnetic impurities belong to. 




FIG. 1: Range functions in the presence of spin-orbit interac- 
tion. The curve labelled with 1 represents \w** (R)\, whereas 
the curve labelled with 2 corresponds to \w"°(R)\. 



As follows from Eq. (10), the gapless fermions (A = 0) 
mediate the indirect exchange interaction that is also 3D- 
isotropic, but with the power-law range function, 



<*W = "MM* {F ° Sss ' +Fl(1 - 6ss ' )} ' (16) 



where 



F = / x*K$(x)dx « 0.31, 



Fi= x 2 K 2 (x)dx ps 0.93. 
Jo 



(17) 



The interaction given by Eq. (16) is the term similar to 
that derived in Ref. 0. It should be noted, however, 
that the gap due to spin-orbit interaction induces an ex- 
ponential decay of the range function, even in samples 
with low density of non-magnetic impurities. 

The indirect exchange interaction is ferromagnetic. 
The numerical results in Fig. 1 show the absolute value of 
the interaction energy when the interacting moments are 
in the same or different sublattices. In our estimations we 
used the carbon-carbon distance a — 1.42 A, the Fermi 
velocity hv = 5.7 eV A, g ~ 1 eV$& and A ~ 1.2 K. 13 
From the unit cell area of S = 3a 2 -\/3/2 containing two 
carbon atoms, one can estimate the coupling coefficient 
go/n~ 2.62 eV A 2 . In order to illustrate the range func- 
tion dependence on R/a, we took into account the scaling 
coefficient a/R = 2.5 x 10~ 5 . 

The numerical results explicitly show the exponential 
decay of the range function with the distance between two 
magnetic impurities. Moreover, they also demonstrate 
that the coupling is stronger when the magnetic moments 
are localized on different sublattices. 
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III. MAGNETIC POLARIZATION DUE TO A 
SINGLE IMPURITY 

The exchange interaction between a magnetic impurity 
and 2D electrons gives rise to a spin polarization of the 
latter ones. The corresponding spin density distribution 
M(i?) can be calculated as a response of the electron 
system to the perturbation due to a single magnetic im- 
purity So at R = 0. As a result, one finds M(i?) in the 
form 

M(R) = -ig Tr J |<rG £ (R)(l + (-l) T r 2 ) 

x (rr-So) G e (-R). (18) 

The assumed exchange interaction between a magnetic 
impurity and 2D electrons is local, see Eq. (3). As a 
result, the spin polarization at a given point of the lattice, 
mediated via the 2D electrons, is sensitive to the location 
of the magnetic impurity (in one sublattice or the other), 
see Eq. (16). 

Performing calculations similar to those presented 
above one finds 

M a {R)= W ^"' SoP . (19) 
5o 

The above result shows that spin polarization has no os- 
cillatory component. Apart from this, one can conclude 
from Eqs. (10) and (19), that the z-component of mag- 
netic density M(i?) is nonzero only in the case when the 
magnetic impurity has a non-vanishing spin component 
perpendicular to the graphene plane. 

IV. EFFECT OF INTERVALLEY TRANSITIONS 

As we have already mentioned above, graphene has 
two nonequivalent valleys in the Brillouin zone, which 
are associated with two Dirac points. Moreover, the ex- 
change coupling between localized moment and 2D band 
electrons assumed in Eq. 3 does not allow for interval- 
ley transitions. However, such transitions are allowed by 
symmetry of the system and should be taken into ac- 
count. To do this we consider an additional term in the 
interaction between a localized moment and band elec- 
trons. This new term may be written in the form 

H *» = — J2 C 1 + M)^*) P» ( S * Ri) , (20) 

n ' 

i 

where g v is the corresponding coupling parameter. 

The above form of interaction between 2D band elec- 
trons and localized moments gives rise to additional 
terms in the indirect exchange interaction. These terms 
are of oscillatory type, ~ cos (2k • R). However, we be- 
lieve that the coupling coefficient g v is much smaller than 
go, so the main results and conclusions described above 
will not be changed by including the contribution follow- 
ing from the interaction (20). 



V. SUMMARY AND CONCLUSIONS 

In this paper we have calculated the indirect exchange 
interaction between two magnetic impurities localized in 
a graphene plane. The interaction is mediated by two- 
dimensional electrons and holes, and acquires the form of 
an isotropic 3D Heisenberg coupling, with the coupling 
parameter being of ferromagnetic type. The exchange 
parameter has an exponentially decaying long distance 
tail determined by the spin-orbit interaction. 

Thus, graphene doped with magnetic impurities may 
become ferromagnetic at a certain impurity concen- 
tration and below a certain temperature. Magnetic 
graphene could be then an ideal system for experimental 
investigations of 2D transport of spin polarized electrons. 
Moreover, one may also expect that magnetic graphene 
could be very useful in spintronics devices, where electron 
spin plays a role comparable to its charged 

We have calculated the exchange interaction for T = 0. 
Strictly speaking, our results refer to T <C A. The effect 
of temperature on the exchange interaction in narrow- 
gap semiconductors has been studied by Rusin^S who 
showed that for the chemical potential in the gap, the 
temperature corrections do not change the main result: 
even for T ^> A, the interaction at large distances is 
exponentially decaying. 

When calculating the exchange interaction we assumed 
the chemical potential fi = 0. The results do not change 
for a finite \i as long as the chemical potential is located 
within the energy gap, which corresponds to the undoped 
material. In principle, the chemical potential in graphene 
can be varied either by a gate voltage or by doping with 
impurities. As a result, the chemical potential can be 
shifted into the conduction or valence bands, and the in- 
teraction between magnetic impurities includes then the 
usual RKKY exchange term via free carriers. As in con- 
ventional metals or doped semiconductors, this gives an 
additional oscillating term with the oscillation period de- 
termined by the Fermi momentum of electrons (holes). 
However, the mechanism of exchange interaction consid- 
ered in this paper is still working for > A. The result- 
ing interaction is then due to both the virtual transitions 
of electrons through the gap and the excitations of real 
electron-hole pairs in the vicinity of the Fermi energy. 

Magnitude of SO interaction in graphene was re- 
cently a subject of extensive discussion in the relevant 
Iiteratureii2i22i2i A similar problem was also discussed 
in the case of carbon nanotubesi22i2&2i Some results in- 
dicate that the SO interaction in graphene can be very 
small. However, the key point is that this interaction 
is nonzero, which results in a finite exchange interaction 
length. As we have shown in this paper, the nonzero 
SO interaction leads to some anisotropy of the magnetic 
interaction between magnetic ions. 

The energy gap in graphene can be also induced by 
excitonic effects thus being not related to the spin 
degrees of freedom. As a result, the exchange interaction 
mediated by virtual transitions through the excitonic gap 
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A exc acquires the finite length hv/A exc . The magnetic 
interaction is then strictly isotropic. 
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